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Data from the BRAHMS collaboration [jl|] show that antiparticle to particle ratios have a max- 
imum at mid-rapidity and slowly decrease as the rapidity becomes larger. It was emphasized very 
eloquently by Rohrich ^ that the particle ratios at large rapidities are consistent with those mea- 
sured at the SPS energies. This opens the possibility to compare measurements for e.g. the K + /n + 
ratio at high rapidities and check them with the corresponding values measured in the energy scan 
at the SPS, thus complementing information about the rapid variation of this ratio as a function of 
beam energy. The sharp variation in this ratio with beam energy remains a mystery for most models 
and might indicate interesting dynamics possibly even linked to the appearance of the critical point 
at intermediate values of the baryon chemical potential in lattice quantum chromodynamics [J3|, Q]. 

A change in particle - anti-particle ratios with longitudinal momentum has been observed in 
many collisions, it was first discussed two decades ago in an analysis of p — p and p — N data [Qj. 

The analysis of particle multiplicities in heavy ion collisions has shown overwhelming evi- 
dence for chemical equilibrium in the final state except for particles carrying strangeness which are 
suppressed at lower beam energies; however, their relative yields fulfill statistical equilibrium. A 
summary as of 2007, combining the results from many different groups, is shown in Fig. |TJ Except 
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Figure 1: Temperature vs. as determined from heavy ion collisions at different beam energies. The lower 
AGS points are based on a preliminary analysis of An data. The RHIC point at 62.4 GeV was obtained by J. 
Takahashi [g]. See [Q, §] for more details. 

for particle multiplicities at RHIC energies, all data in Fig. [j] use integrated particle yields, the very 
systematic change of thermal parameters over the full range of beam energies is one of the most 
impressive features of relativistic ion collisions to date. It is now possible to use the thermal model 
to make reliable predictions for particle multiplicities at LHC energies []|] and to determine which 
beam energy will lead to the highest baryon density at freeze-out [|i0|]. With chemical equilibrium 
thus firmly established, we focus on other properties, in particular, since the rapidity distributions 
of identified particles is now becoming available also at RHIC energies [Q], it is now possible to 
determine the rapidity dependence of thermal parameters. A first analysis was done by Stiles and 
Murray [JTTJ] for the data obtained by the BRAHMS collaboration at 200 GeV [g]. This shows a 
clear dependence of the baryon chemical potential on rapidity due to the changing p/p ratio. A 



thorough analysis of the rapidity dependence was recently done in Ref. [12, O] using a model 
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based on a single freeze-out temperature. A first report of our results was presented elsewhere [|14|] . 

The general procedure is as follows: the rapidity axis is populated with fireballs following a 
gaussian distribution function given by p(yFE>) where ypB is the rapidity of the fireball. 

Particles will appear when the fireball freezes out and will follow a thermal distribution cen- 
tered around the position of the fireball The momentum distribution of hadron i is then calculated 
from the distribution of fireballs as given by Eq. [jl|] along the rapidity axis as follows 

d 3 Ni f°° d 3 N', 
Ei ^P = L P (y™)E i7 ^(y-y FB ) dy FB (2) 

where Ej^J^ is the distribution of hadrons from a single fireball. The temperature T and the baryon 
chemical potential /j.b will depend on the rapidity of the fireball and are not assumed to be constant. 

We have included resonance decays in the final distribution and assumed they decay isotropi- 
cally. 

An important parameter is the width of the distribution. For the RHIC data at 200 GeV this 
was determined from the 7T + 's as these are very sensitive to the value of a and less to variations 
in the baryon chemical potential. The width of the distribution a = 2.183 is compatible with the 
values quoted by the BRAHMS collaboration [@], e.g. c n + = 2.25 ±0.02 and a K = 2.29 ±0.02. 
This is shown in Fig. The hadrons described by Eq. [§] are mainly hadronic resonances. Only a 
fraction of these are stable under strong interactions and most of them decay into stable hadrons at 
chemical freeze-out, hence the need to implement multi-particle decays. 

We assume that the temperature T and the chemical potential /ig are always related via the 
freeze-out curve as given in Fig. §\§ ; if the temperature varies along the rapidity axis, then also the 
chemical potential will vary. Thus a decrease in the temperature of the fireball will be accompanied 
by an increase in the baryon chemical potential. In other words, we assume a universality of 
the chemical freeze-out condition. This relationship between temperature and baryon chemical 
potential is very reasonable since all particle abundances measured so far follow it. Once the width 
of the distribution of fireballs has been fixed, we can go on with the dependence of the baryon 
chemical potential on the rapidity of the fireball (units are in GeV) at RHIC: 

Hb = 0.025 + 0.01 lyp B (3) 

This is shown graphically in Fig. [g]. 

A comparison of the resulting net proton distribution p — p with RHIC data at 200 GeV is 
shown in Fig. [Qj. 

The variation of the temperature along the rapidity axis is shown in Fig. The temperature 
is maximal at mid-rapidity and gradually decreases towards higher (absolute) values of the rapidity. 
Note that the temperature varies less than 2 MeV over the rapidity interval. 

The situation at the highest SPS energy is more difficult to describe because the changes of 
the thermal variables with rapidity are much larger. First of all we need the distribution of fireballs 
at SPS energies. This is determined by the variable a appearing in Eq. ([!]) which we fixed using 
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Figure 2: Fit to the pion distribution as measured by the BRAHMS collaboration. 
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Figure 3: Values of the baryon chemical potential as a function of rapidity at the highest RHIC energy. 

the distribution of pions. The results obtained at various beam energies are shown in Fig. [jg]. The 
baryon chemical potential can approximately be described by: 

Hb = 0.237 + 0.01 1 yj? B . (4) 

This is shown graphically in Fig. [JsJ] . The corresponding change in the temperature as a function of 
rapidity is shown in Fig. Whereas there is almost no change in the temperature in the rapidity 
interval under consideration (the change is less than 2 MeV at RHIC), at the highest SPS energy 
the change in temperature is much larger and ranges from about 160 MeV down to about 120 MeV. 
The comparison between the resulting net proton distribution p — p with SPS data at 17.2 GeV is 
shown in Fig. [|7|]. 

In summary, particle yields measured in heavy ion collisions show an overwhelming evidence 
for chemical equilibrium at all beam energies. The rapidity dependence of the thermal parameters 
T and /ig can now be determined over a wide range of rapidities and show a systematic behavior 
towards an increase in as the rapidity is increased. 
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Figure 4: The p — p distribution as a function of rapidity. 
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Figure 5: Values of the chemical freeze-out temperature function of rapidity at the highest RHIC energy. 
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Figure 6: Width of the Gaussian distribution, a, used to describe the distribution of fireballs for various 
beam energies. 
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Figure 7: The p — p distribution as a function of rapidity at SPS. 
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Figure 8: Values of the baryon chemical potential as a function of rapidity at the highest SPS energy. 
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Figure 9: Values of the chemical freeze-out temperature as a function of rapidity at the highest SPS energy. 
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